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Abstract 

We examine the finite volume spectrum and boundary energy in boundary sine- 
Gordon theory, based on our recent results obtained by closing the boundary boot- 
strap. The spectrum and the reflection factors are checked against truncated confor- 
mal space, together with a (still unpublished) prediction by Al.B. Zamolodchikov for 
the boundary energy and the relation between the parameters of the scattering ampli- 
tudes and of the perturbed CFT Hamiltonian. In addition, a derivation of Zamolod- 
chikov's formulae is given. We find an entirely consistent picture and strong evidence 
for the validity of the conjectured spectrum and scattering amplitudes, which together 
give a complete description of the boundary sine-Gordon theory on mass shell. 

Introduction 



Sine-Gordon field theory is one of the most important quantum field theoretic models 
with numerous applications ranging from particle theoretic problems to condensed matter 
systems, and one which has played a central role in our understanding of 1 -|- 1 dimensional 
field theories. A crucial property of the model is integrability, which permits an exact 
analytic determination of many of its physical properties and characteristic quantities. 
Integrability can also be preserved in the presence of a boundary if suitable boundary 
conditions are imposed 

In this paper, continuing our work started in ||2|, 0|, we investigate sine-Gordon field 
theory on the half-line and on a finite volume interval, with integrable boundary condi- 
tions. It was first pointed out by Ghoshal and Zamolodchikov [|] that the most general 
integrable boundary potential depends on two parameters. They also introduced the notion 
of 'boundary crossing unitarity', and combining it with the boundary version of the Yang 
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Baxter equations they were able to determine soliton reflection factors on the boundary; 
later Ghoshal completed this work by determining the breather reflection factors using 
a boundary bootstrap equation first proposed by Fring and Koberle 0. 

The results of Ghoshal and Zamolodchikov concerned only the reflection factors on 
the ground state boundary, although they already noticed that there are poles in the 
amplitudes which signal the existence of excited boundary states. The first (partial) results 
on the spectrum of these boundary states were obtained by Saleur and Skorik for Dirichlet 
boundary conditions 0. However, they did not take into account the boundary analogue 
of the Coleman-Thun mechanism, the importance of which was first emphasized by Dorey 
et. al. 0. Using this mechanism Mattsson and Dorey were able to close the bootstrap in 
the Dirichlet case and determine the complete spectrum and the reflection factors on the 
excited boundary states 0. Recently we used their ideas to obtain the spectrum of excited 
boundary states and their reflection factors for the Neumann boundary condition |3 and 
then for the general two-parameter family of integrable boundary conditions 0]. For the 
Neumann case, extensive checks were performed using a boundary version of the so-called 
Truncated Conformal Space Approach (TCSA) ||10|, |1T|; for the generic case, however, these 
checks were not carried out at that time. 

Another interesting problem is that of the boundary energy. Namely, the boundary 
contributes a volume independent (constant) term to the free energy, in addition to the 
bulk energy density which gives a term proportional to the spatial volume. Just as in the 
case of the bulk energy density, the boundary energy in general QFT is not a universal 
quantity. However, in perturbed conformal fleld theories there is a preferred normalization^ 
of the Hamiltonian which gives a unique deflnition for both the bulk and the boundary 
contributions. Therefore, this boundary energy is an interesting quantity to compute. For 
Dirichlet boundary condition it was obtained by Leclair et al. in ||T^. A few years ago Al. 



B. Zamolodchikov presented a result for general integrable boundary conditions ||T^ . 

One crucial ingredient, that is needed e.g. for a TCSA check of the spectrum and 
reflection factors for the general integrable boundary conditions, is a relation between the 
ultraviolet (UV) parameters that appear in the perturbed CFT Hamiltonian and the in- 
frared (IR) parameters in the reflection factors. This relation was also obtained by Al. 



B. Zamolodchikov ||13|. Using his result, we perform an extensive check of the spectrum, 
boundary energy and reflection factors of boundary sine-Gordon theory. This provides 
strong evidence that all the results mentioned above form a consistent and complete de- 
scription of the boundary sine-Gordon theory on mass shell (i.e. spectrum and scattering 
amplitudes). 

The paper is organized as follows. In Section | we recall the results on the boundary 
bootstrap in boundary sine-Gordon theory. Section ^ describes Zamolodchikov's formulae 
on the UV-IR relation and the boundary energy. These formulae were presented at some 
seminars, but have not been published; however, we could get some notes taken by the 

^In this preferred normalization, the perturbing bulk and boundary operators transform homogeneously 
under scale transformations. 
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audience.0 In these notes we found several misprints; in order to determine the correct form 
of the formulae (which is of utmost importance in order to make the comparison to TCSA), 
we rederive here the boundary energy using the thermodynamic Bethe Ansatz (TBA) and 
then check the UV-IR relation using the exact vacuum expectation values of boundary 
fields conjectured by Fateev, Zamolodchikov and Zamolodchikov 10 In Section ^ we 
describe the results coming from TCSA for generic (non Dirichlet) boundary conditions, 
while in Section |] we present the results for Dirichlet boundary conditions, which are a 
singular limit of the generic case and so the TCSA must be set up differently. The paper 
ends with some brief conclusions and an outlook in Section ||. 



2 Boundary bootstrap in sine-Gordon theory 

Boundary sine-Gordon theory is defined by the action 

/oo / 
dt i dx 
CO \J —QO 

where ^{x, t) is a real scalar field and Mq, 0o are the two parameters characterizing the 
boundary condition: 

a.$(x,t)U=o = -Mo^sin (^^($(0,t) -0o)) • (2.2) 

Ghoshal and Zamolodchikov showed that the above model is integrable |||] and that the 
above boundary potential is the most general that permits the existence of higher spin 
conserved charges. 
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+ Mocos^($(0, t)-0o; 



(2.1) 



2.1 Bulk scattering properties 

In the bulk sine-Gordon model the particle spectrum consists of the soliton s, the antisoliton 
s , and the breathers -B" which appear as bound states in the ss scattering amplitude S^l . 
As a consequence of the integrable nature of the model any scattering amplitude factorizes 
into a product of two particle scattering amplitudes, from which the independent ones in 
the purely solitonic sector are 



a{u) = S^^{u) = S (u) 

b{u) = Slz{u) = 51+ (n) 



n 

1=1 



Alt ^ 



T{2{1 - 1)X - f)T{2lX + 1 - f 



r((2/-l)A-^)r((2/-l)A+l-^ 



-/{u -u) 



sin(Au) 



ciu 



sin(A(7r — u)) 

sin(A7r) 
sin(A(7r — u)) 



a{u) 



5 



(2.3) 



^We thank P. Dorey and G. Watts for making the notes available to us. 

^The derivation presented here is very similar to the way ALB. Zamolodchikov arrived to the formulae 



( P ■ 2 - 3^ ) in Section 3.1, according to the hints he gave in his seminars. 
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where the parameter A is determined by the sine-Gordon coupling constant 



A-^^ 1 



(2.4) 



and u = —iO denotes the purely imaginary rapidity. The other scattering amplitudes can 
be described in terms of the functions 



{y} 



W = ■ )l li , {y}{-y} = 1 , {y + 2A} = {-y} 



sm 



f U XTT \ 

\2 2 ) 



as follows. For the scattering of the breathers -B" and -B™ with n > m and relative rapidity 
u the amplitude takes the form 

^"'"(m) = Slil^iiu) = {n + m- l}{n + m - 3} . . . {n - m + 3}{n - m + 1} , 

while for the scattering of the soliton (antisoliton) and 5" 

{1 + A} if is even 



S^{u) = S+^H = SZlliu) = {n-l + X}{n - 3 + A} 



2.2 Ground state reflection factors 



if n is odd 



The most general reflection factor - modulo CDD-type factors - of the soliton antisoliton 
multiplet \s,s) on the ground state boundary, denoted by |), satisfying the boundary 
versions of the Yang Baxter, unitarity and crossing equations was found by Ghoshal and 
Zamolodchikov H: 



i?(?7, -d, u) 



P+{ri,i3,u) Qiri,i3,u) 

Q{rj,{},u) P^{rj,{},u) 

P^{v,^,u) Qo{u) 

Qoiu) Pq {ri,^,u) 



Ro{u) 



a{r], u) a{i^, u) 



Qoiu) 



cos{ri) cosh('i9) 
cos(Am) cos(?7) cosh('i9) =i= sin(A-u) sin(?7) sinh('(9) 
— sin(A'u) cos(A'u) 



(2.5) 



where rj and -(9 are the two real parameters characterizing the solution 

r(4/A - ^)r(4A(/ - 1) + 1 - 



Ro{u) = n 



1=1 



r((4/ - 3)A - ^)r((4/ - 1)A + 1 - 2^ 
is the boundary condition independent part and 

r(| + f + (2/ - i)A - ^)r(i - £ + (2/ - i)A - ^ 



ct(x, u) 



cosx 



cos(x + Am) 



n 

1=1 



r(i - f + {21 - 2)A - ^)r(i + I + 2/A - ^) 



-/{u -u) 
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describes the boundary condition dependence. Note that the topological charge may be 
changed by two in these reflections, thus the parity of the soliton number is conserved. 

As a consequence of the bootstrap equations the breather reflection factors share 
the structure of the solitonic ones, [|]: 

R^''\ri,^,u) = R^^\u)S^''\ri,u)S^''\i'&,u) , (2.6) 

where 

\2\^ 2) 1=1 ilA + 2 J «=0 Utt 2 2A ) 

(n) 

In general Rq describes the boundary independent properties and the other factors give 
the boundary dependent ones. 

2.3 The spectrum of boundary bound states and the associated 
reflection factors 

In the general case, the spectrum of boundary bound states was derived in 0. It is a 
straightforward generalization of the spectrum in the Dirichlet limit previously obtained by 
Mattsson and Dorey |||]. The states can be labeled by a sequence of integers \ni, n2, . . . , n^). 
Such a state exists whenever the 



TT 

- > Z^ni > > Z/„3 > > . . . > 



condition holds, where 



n (2n + l)7r , n (2n - l)7r 

Un = ; and W„ = TT — ; 

A 2A A 2A 

denote the location of certain poles in a{rj,u). The mass of such a state (i.e. its energy 
above the ground state) is 

"^|ni,n2,...,nfe> = M ^ COs(z/„J + M ^ COs(w„J , (2.8) 
i odd i even 

where M is the soliton mass. The reflection factors of the various particles on these 
boundary states depend on whether k is even or odd. When k is even, the reflection 
factors take the form 

Q|ni,n2,...,nfe>(r?,^?,M) = <5(r?,^?,M) JJ^^^^^'") 11 ' 

i odd i even 

and 

^|ni,n2,...,nfe>(^''^'^) = ^^(^'^'^) H ""«(^'^) H ' 

i odd i even 
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for the solitonic processes, where 

n 

«n(r7,M) = n{2(^-/)} ; f] = Ti{\ + l)-r, ■ 
1=1 

For the breather reflection factors the analogous formula is 

<tn„...,n,)(^'^'«) = R^''\v,^,n) n hlXv.u) n hl^V.u) (2.9) 

i odd i even 

where now 

hl{'q,u)= W \^-\ + n-2l\\^ + \-n~2{k+l-l)\ . (2.10) 

In the case when k is odd, the same formulae apply if in the P^, Q and ground state 
reflection factors the rj ^ f] and s ^ s changes are made. 

3 Boundary energy and UV-IR relation in sine- Gordon 
theory 

3.1 Zamolodchikov's formulae 

As mentioned in the introduction, recently Al. B. Zamolodchikov presented (but not yet 
published) ||13| a formula for the relation between the UV and the IR parameters in the 
sine-Gordon model. To set the conventions for this relation, consider boundary sine-Gordon 
theory as a joint bulk and boundary perturbation of the c = 1 free boson with Neumann 
boundary conditions (perturbed conformal field theory, pCFT): 

dt dx : cos/5$(x, t) : +jx j dt : cos W<l>(0, t) - 0o) : (3.1) 

-oo J —oo J —oo ^ 

where the colons denote the standard CFT normal ordering, which defines the normaliza- 
tion of the operators and of the coupling constants. The couplings /x and jl have nontrivial 
dimensions; 



[/i] = [mass]^-^'*'', [p\ = [mass]'-''^ h 



see the section on TCSA for more details. With these conventions the UV-IR relation Q is 

cos cosh = - — cos 



TT / \87r J yUcrit V 2 



sin(— Jsinh(^— J = -^sin( — ) , (3.2) 



''A similar relation was derived by Corrigan and Taormina ||T^ for sinh-Gordon theory, however, their 
normalization of the coupling constants is different from the one natural in the perturbed CFT framework. 
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where 



/^crit 



2/i 



sm 



/3= 



(3.3) 



Zamolodchikov also gave the boundary energy as 
M 



^/ nx M f fr]\ , 1 / TT \ 

Eiri, v) = cos — + cosh — cos — 1 

' 2cos|^ V 2 V2A 



1 / TT 

- sm — 

2 V2A 



(3.4) 



3.2 Derivation of the boundary energy from TBA 



In an integrable boundary theory with one scalar particle of mass m only, one can write 
down the TBA equation for the ground state energy on a strip with spatial volume L and 
integrable boundary conditions a and b at the two ends. The equation is of the form p^ : 

/oo in/ . , 

— ^ (^ - 0') log (l + e"'^''^) (3.5) 

where / = mL is the dimensionless volume parameter. The kernel is expressed in terms of 
the two-body S'-matrix S{ff) as 

d 



while 



Kh{Q) = - log 



Rn ( ^1 Rh i 
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where Ra (9) and Rb (9) are the reflection factors for the two ends. From the solution e (9) 
of the TBA equation the ground state energy can be calculated using the formula 



E{L) = E],u\kL + E] 



boundary 



TTc{mL) 
2AL 



c{l) = — / d9L{9) cosh9 , (3.6) 

J-oc 



where L{9) is the usual short hand notation L{9) = log (l + e"^'-^-*) . It is well-known 
that no such TBA equation (or, for that matter, a flnite system of TBA equations) can be 
written for sine-Gordon theory as a result of the nondiagonal bulk and boundary scattering 
of the solitons (except for special values of the parameters). Therefore, our approach is to 
calculate the boundary energy for sinh-Gordon theory and then analytically continue back 
to the sine-Gordon case. This is known to work e.g. for 5'-matrices, form factors and many 
other quantities, and so we simply assume it works for the boundary energy as well. 
Consider therefore the boundary energy in boundary sinh-Gordon theory 



A. 



hG 



dt 



dx 



1 TtX 

-d^<^d^'<l> - -^cosh&<l> 



MoCosh-($-0o; 
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which can be considered as the analytic continuation of the boundary sine-Gordon model 
by substituting h = i(3 (and changing the convention for the sign of Mq). Then from 

and, as a result, A is negative for the sinh-Gordon case. Note that the analytic continuation 
is through the point \ = oo (complex infinity), therefore for the purposes of relating 
physical quantities between the two models the natural variable is A~^. 

We now proceed to the calculation of the boundary energy. A similar calculation was 



performed by Dorey at al. for the scaling Lee- Yang case. They presented the general 
idea with enough hints to reconstruct the method, but for the sake of completeness we 
write down the details for the interested reader. It is based on Zamolodchikov's method 



for obtaining the bulk energy from the TBA with periodic boundary conditions |[T^ . 

Suppose for simplicity that the boundary conditions a and h are identical and so = kaa 
is even. Then in general the functions k and ip have the following asymptotic behaviour 

k{e) ~ /to + Ae-i^'i + . . . 

^{d) ~ Ce-I^l + ... (3.7) 

for \6\ oo, where k^, A and C are real constants. 
The 'kink' functions, defined as 

e±{e) = lime (^±logy 

satisfy the 'kink' equation 

(9) = e^' + ko- 1^ {9 - 9') log (l + e-±(^')) 

and are related as 

e (9) = e+ i-9) . 
Let us also introduce the following definitions 

L±(^) =log(l + e-^±('^)) 

and define the asymptotic values 

eo = e+{-oo) , Lo = L+(-oo) = log (l + e"^°) 

which satisfy the standard 'plateau' equation 

/CO in 
^viO) . (3. 
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Our aim is to expand c(/) around / = 0. To calculate the first few terms, it is convenient 
to define the functions 6 and L in the following way: 

m = L+(e-\ogj^+L+(^-9-\ogj^+L{9)-Lo. (3.9) 

They satisfy 

/oo jnl 
^(9 -9') 1(9') , (3.10) 

6{9) , L{9) ^ asl 0. 

We can then rewrite 

c{l) = ^ d9e'L+{9) + ^ d9e-'^ + ^ I d9L{9)cosh9 



The first term gives the UV central charge and can be calculated using the standard 
dilogarithm sum rules. The second term is the (anti) bulk energy density, that can be 
calculated self-consistently by examining the 9 ^ —oo asymptotics of the integrand [pT^j: 

dL+ 1 de+ 1 Oea 



89 l + e=+(^) 89 1 + 89 



for 9 — s> — oo 



The terms proportional to must cancel for the integral to converge on its lower bound. 
Using the kink equation and the asymptotics of y?, to leading order 



89 \ 2n .L^ 89' 
from which 

-W=c- 

In the perturbed conformal field theory formalism, the ground state energy can be expanded 
as 

oo 

TT 



n=0 

SO the terms linear in L must cancel from ( ^.61) . Therefore we obtain the bulk energy 
density as 

The third term can be rewritten using that L{9) = L {—9): 



/oo /"OO 
d9 L (9) cosh 9= d9L{9) e 
■oo J — oo 
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After a partially integration, it can be seen that once again, the integral is convergent 
if terms proportional to cancel in Using equations (|3.9| ) this is equivalent to can- 
cellation of all terms proportional in 6, at least to leading order in /. From (|3.10| ) we 
obtain 

' dO'L {6" 

— oo 



6{e)=e^ [A- — I dO'L [9') e 



from which we obtain (to leading order) 

2tiA 



deL{e) e- 



c 



None of the subleading terms contains any contribution which are independent of the 
volume and therefore in ( |3.6| ) -^boundary must cancel against this particular term, leading to 

A 



-^boundary ~ • 



2C 



3.3 The sinh-Gordon case 



In sinh-Gordon theory the two particle S'-matrix can be written as (remember, that in our 
convention A is negative in its physical range): 

, ^ sinh 9 + isin^ , ^ 

S (0) = -TTT--^ ■ (3-11) 



sinh 9 — i sin j 



As a result, the TBA kernel is 



, , 2 cosh 6^ sin Y vr ia\ i i\a\\ 

m = - . ■ 2\ - -4sin-e-l^l+0 e-2|^l 

smh 9 + sm j A 

and so we get 

(7 = -4 sm - . 

A 

The integral N takes the value 

2n^^ ' \ -1 for 3?eA > 



which means that the plateau equation ( p. 8]) has the solution 



e-.o _ ^ 



1 + e ^0 sign 3f?e A 



Note that there is no real solution for A < 0, fco < 0. This peculiarity of the sinh-Gordon 
TBA equation was already noted by Al.B. Zamolodchikov in the case of periodic boundary 
condition ||18|. We simply assume that we are working for parameter values for which such 
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a solution exists, so the considerations of the previous subsection apply. This is always the 
case for 3?e A > 0, which is however not a physical range of the parameter A in sinh-Gordon 
theory. Therefore we treat the sinh-Gordon TBA in this range as a mathematical problem 
only, without a corresponding physical field theory (except for the case A = 3/2, see later). 
We further assume that all physical quantities that we wish to calculate are meromorphic 
functions of A~^ and so they have a unique analytic continuation to the values of A~^ that 
we are interested in0. Note that this argument is not a proper derivation; however, for 
the time being this is the only way we can arrive at the desired result, and we show that 
the results fit with the bootstrap spectrum, TCSA data and known results from previous 
literature. 

E.g. for the bulk energy density we obtain 

2 

pshG _ ^ 
-^bulk — -—■ - ■ 



This is meromorphic in A ^ and so we trust that it is the true bulk energy constant of the 



sinh-Gordon theory in the regime A < 0. Furthermore, it is equal to the known result ||19 
Now we can try and continue this result to the sine-Gordon regime A > 0. Under this 
continuation the sinh-Gordon particle is identified with the first breather of sine-Gordon 
theory and so we have 

71 

m = 2Msin — , 
2A ' 

where M is the soliton mass. We then obtain 

C _ M2 TT 
-^bulk — ^ ^ ) 

which is the correct bulk energy density of sine-Gordon theory ||2^. Hence the above 
method of continuation works for the bulk energy constant. 

Now let us calculate the boundary energy. From eqns. (|2.6| , |2.7|) , the refiection factor 
of the first breather in sine-Gordon theory can be written as 

where r] and d parametrize the boundary conditions. The sinh-Gordon refiection factor can 
be obtained by continuing the refiection factor to negative values of A~^ (for sinh-Gordon 
theory, rj is real and d is purely imaginary, while for sine-Gordon theory both parameters 
are real). Putting the same boundary condition on the two boundaries of the strip (with 
the same values of and rf) we obtain 



^It is clear that the relevant variable to consider is A ^ because the continuation in the coupling goes 
through the value /? = which corresponds to A = c» 
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where E{ri, {}) is the energy of a single boundary. The term k {6) in the TBA equation 
(PD is 

k{e) = - log [K {6) K {-6)] , K{9) = (z| - . 
Using the identity 

cosh 9 + sin nx 
* 2 * 2 cosh 9 — sin vrx ' 

we get 

-log[(x)i|+e(x)if-e] ~ -4sin7rxe-l^l + 0(e'2|e|^ ^ 

Note that and A in ( p.7| ) can be calculated additively from the asymptotics of the 
contribution of a single (x) block above. As a result, ko = and so the plateau eqn. 
( p.8| ) has no solution in the sinh-Gordon regime A < 0, thus the analytic continuation 
described above cannot be avoided. Putting the ingredients together, the boundary energy 
in sinh-Gordon theory takes the form 

It is now easy to recover Zamolodchikov's formula (|3^) for the boundary energy in sine- 
Gordon theory. 

As an immediate check on this calculation, we wish to note that for A = | the ^-matrix 
( p.llj ) is identical to that of the scaling Lee- Yang model, and the reflection factors of the 
scaling Lee- Yang model corresponding to integrable boundary conditions are reproduced 
by specifying some complex values for rj and {}. It can be easily checked that the formula 
( 3.13| ) reproduces correctly the results of Dorey et al. ||TT| . 

3.4 Special cases 

Since we obtained the boundary energy of sine-Gordon/sinh-Gordon theory under some 
non trivial assumptions we check the results in some known cases. 

3.4.1 Dirichlet boundary conditions 

Dirichlet boundary conditions correspond to the limit /i — > cxd in (^), which leads to 

2n 

$(0, t) = 00 mod — . 



The reflection factor of the first breather can be obtained as the 'i? — > oo limit of ( p. 12 



(1) (±4-1) f^-i) 
\2X~ 2/e \nX~ 2/9 
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The derivation of the previous subsection then gives the boundary energy 

EdM = I cos (^) - - cos (^) + - sin ( ^] - - ) , (3.14) 

^" 2cos|^ V 2 V2A/ 2 V2A/ 2/' ^ ' 



which is exactly identical to the formula obtained by Leclair et al. in The parameter 
rj is related to 0o in the following way 

r] = 7c{X + l) — , 
Zn 

which was conjectured by Ghoshal and Zamolodchikov 0, and is a straightforward conse- 
quence of eqns. ( |3.2| ) as well. 

Note that Ed^t]) cannot be obtained as the {} ^ oo limit of the general boundary 
energy eqn. (|3.4|) . The reason is clear: the boundary potential is normalized in different 
ways in the two cases: classically to obtain finite energy in the Dirichlet limit one has to 
add Mo to the general — Mq cos (f ($ — (po)) boundary potential. Clearly in the quantum 
case, when the boundary vertex operator has a non trivial dimension, we can not simply 
subtract p, from Eijj.'d). Since the quantity we subtract must have the dimension of mass 
and should depend on /2, it must be proportional to = The question 

is whether we can make this subtraction such that in the ^ oo limit the leading term 
cancels and the constant terms just reproduce Eniji). The UV-IR relations, eqn. ( fj.2[ - P7^ ) 
guarantee, that 



Thus, upon using the bulk mass gap relation (cf. eqn. ([4.2|)), becomes proportional 



to Me^f^ up to exponentially small terms for •§ ^ oo. Therefore, if we subtract this term 
with an appropriate coefficient then in the Dirichlet limit the surviving constant terms 
exactly reproduce (|3.14|) . 



3.4.2 The first excited state 

It was noted in [^] (for Dirichlet boundary condition) and in (for the general case) that 
continuing analytically 

7] 7r(A + 1) - ?7 

the role of the boundary ground state |) and the boundary first excited state |0) are 
interchanged. Therefore we can calculate the energy difference between these two states 
from the formula for the boundary energy, eqn. (|3.4|) . The result is 

(7r(A + 1) - 77, ^?) - i?(r/, ^)=M cos (| - ^) 

which exactly equals the prediction of the bootstrap, i.e. 

E\Q) — E\) = Mcosz/q 

that follows from eqn. (|2.8|). 
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3.5 UV-IR relations and vacuum expectation values (VEVs) 

As it is well known in the bulk sine-Gordon theory there is a relation among the following 
three exactly calculable quantities: the ground state energy density, the dimensionless 
constant entering the mass gap relation connecting the UV and IR parameters, and the 
VEV of the exponential field (e*^**^^^) [2^. This relation is such that knowing any two 



of these quantities determines the third one. It generalizes to sine-Gordon theory with 
boundaries, where it connects the boundary energy, the UV-IR relations (|3.2[ - p.3|) , and 
the VEV of the boundary field (e^t*(°)) in a similar way. As the VEV of the boundary 
operators has been determined by Fateev, Zamolodchikov and Zamolodchikov (FZZ), we 
can use it to show that the UV-IR relations, (|3.2[ - |3.3|) and the boundary energy, (|3.4|) , are 



indeed consistent with the VEVs given in 14]. For simplicity we consider only the special 



case when 0o = 0, as this case already illustrates the point. (More precisely the condition 
00 = can be satisfied in two different ways 0: either by = or by = 0, and we 
consider the former possibility). 

Writing the functional integral representation of the partition function Zgh = Tie'^^'^''^^^ 
on a cylinder of length R and circumference L with boundary states a and b on the bound- 
ary circles and considering the R ^ oo limit (when Zab ~ e~^^°'''^^^) one readily derives 
that in this limit the ground state energy Eaa satisfies 



djl 



(e^t*(o))^_G'(/5,/,). (3.15) 



(In writing this equation we assumed that G{P, jx) = G{—[3, p,)). Since for {} = the ground 
state energy depends on /2 only through the rj parameter appearing in the boundary energy, 
eqn. ( p.l5| ) actually determines the dependence of 77 on /i. Furthermore, both sides of ( p.l5| ) 
can be integrated to obtain the following expression for the boundary energy 

E{ri) = - j dflG{P,fl). (3.16) 



What we show below is that using the FZZ expression for G{(3, fi) on the r.h.s. gives (^) 
for the boundary energy only if ( |3 .21 - ^.31 ) hold. 

The expression given in ||l^ for G'(/5, p) depends on jl through a parameter 2;, which, 
for 00 = 0, we take to be pure imaginary z = iZ [Z real): 

cos2(7rZ) = |^sin(^^^ . (3.17) 

Then ^ 

( p/ A N \ 2A 
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where go and gs are given by the integral representations^: 



loggo{f3) = I J 



2sinh(t/(A 



gtA/(2A+2) ^^gi^ ( t_ I ^^^^ 



sinh(t) sinh(tA/(A + 1)) 

oo 

dt 2 sinh(t/(A + 1) sinh^iZt) 



t 



2A + 2 



A + 1 



\oggs{P,Z) = - 



t sinh(t)sinh(tA/(A + 1)) 



The integrals appearing here can be computed analytically after some efforts. Finally, 
expressing fi in terms of the soliton mass M via (|4.2| ), and converting the integral over /2 
to an integral over nZ by using ( |3.17|) , after some algebra one finds 



M 



2cos(^) 



cos 



Zn{X + l) 
A 



+ /(A). 



This agrees with the boundary energy, (|3.4|) , if Ztt = i.e. when eqn. (p.l7|) becomes 
identical to (the = case of) (^^^). 



4 TCSA: general integrable boundary condition 
4.1 TCSA for the boundary sine-Gordon model 

First we describe the Hamiltonian of boundary sine-Gordon model (BSG) living on the 
line segment < x < L as that of a bulk and boundary perturbed free boson with suitable 
boundary conditions. This is the starting point of the TCSA analysis. 

The basic idea of TCSA is to describe certain 2d models in finite volume as relevant 
perturbations of their ultraviolet limiting CFT-s [|T^]. If we consider boundary field theo- 
ries, then the CFT-s in the ultraviolet are in fact boundary CFT-s. The use of TCSA to 
investigate boundary theories was advocated in [ll|, [211 . 

As the bulk SG model can be successfully described in TCSA as a perturbation of the 



c = 1 free boson it is natural to expect that the various BSG models are appropriate 
perturbations of c = 1 theories with Neumann or Dirichlet boundary conditions. Therefore 
we take the strip < x < L and consider the following perturbations of the free boson, as 



described in detail in 



oo JO 



S = I I ( -(9^$a^$ + /icos(/5$) ) dxdt + 



j ^/io cos ($5 - 0o)^ + fj'L COS ($B - (Pl)^^ dt . 



^Note that the integral for log 55 contains a factor of 1/2 compared to the expression in [Tij even after 
accounting for the difference between the parameters of this paper and of |Q. Without the inclusion of 
this factor it would be impossible to obtain the correct rj dependence of the boundary energy as in eqn. 
( P.4D . The fact that this factor should be present was later confirmed to us by ALB. Zamolodchikov in a 
private discussion. 
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Here, for finite /I's, Neumann boundary conditions are imposed in the underlying c = 1 
theory on the boundaries, while if any of the jl -s is infinite then the corresponding term is 
absent and the boundary condition in the underlying conformal theory on that boundary 
is Dirichlet. The Hamiltonian of the system can be rewritten in terms of the variables 
associated to the plane using the map {x,it) = ^ ^ z = e*i^: 

^(^)''(e~^^^-vl/^(-l) + e^^*-v[;_^(-l)) . (4.1) 

Here Vp{z,z) = n{z,z) : e*^**-^'^-' : and §_{y) =: e*^**-^'^^ : are the bulk and boundary 
vertex operators and the normal ordering coefficient n{z, z) depends on the boundary 
conditions chosen |0. 

Now the computation of the matrix elements of the bulk and boundary vertex operators 
V±p and \l/±/3/2 (with conformal dimension hf^ = ^) between the vectors of the appropriate 
conform Hilbert spaces is straightforward and the integrals can also be calculated explicitly. 
The TCSA method consists of truncating the Hilbert space at a certain conformal energy 
level i^cut (which is nothing but the eigenvalue of the zeroth Virasoro generator) and 
diagonalizing the Hamiltonian numerically. 

It is important to realize that one has to write separate TCSA programs for checking 
the Dirichlet limit and the general two parameter case. In the Dirichlet case there are no 
relevant operators on the boundary, thus both fiQ and fii must be set to zero, and we can 
have jl-s different from zero only if we perturb a CFT with Neumann boundary condition. 
Therefore we investigate the general two parameter boundary sine-Gordon theory by de- 
scribing it as an appropriately perturbed c = 1 CFT with Neumann boundary conditions 
at both ends. The Hilbert spaces of the c = 1 CFT-s with Dirichlet or Neumann boundary 
conditions at the two ends are rather different: while in the former case it basically consists 
of the vacuum module only, in the latter it is the direct sum of modules built on the highest 
weight vectors carrying the allowed values of the field momentum. 

Let us investigate the general two parameter BSG first. Then the simplest choice 
(i.e. the one resulting in the least complex spectrum which is enough to compare to the 
predictions) is to switch on the boundary perturbation only at one end of the strip. The 
TCSA Hamiltonian for BSG with Neumann boundary condition at one end and perturbed 
Neumann condition, (|2.2|) , at the other, is obtained from (|4.1|) by setting fli = 0, /2o = 
/2 7^ . The spectrum of vertex operators in this case is Vn:{z,z) and \l/™(?/), where r 
is the compactification radius of the free boson of the c = 1 theory in the UV, and n, 

m are integers. These fields are primary under the chiral algebra U{1) (i.e. U{1) affine 
Lie algebra). However the compactification radius must be chosen so that both V±/3 and 
, ^be in the spectrumiFI r = 2\/A7i/P = 2ro. Then V±r are represented as V:2 while \1/ , ^ 

2 r -^2 

^The ^/Att has its origin in the different normaHzations of the SG scalar field <J> and the c = 1 CFT one. 
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as In other words we have to consider the boundary perturbation of the 2-folded 

r 

sine-Gordon model in the sense of 
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We choose our units in terms of the sohton mass M. The bulk coupling fj, is related to 
M by 

fi = K{f3)M'-''^, hp = f^ (4-2) 

where is a dimensionless constant. In the bulk SG, from TBA considerations, the 
exact form of k(/3) was obtained in and we use the same form also here in BSG. Once 



we expressed fi and used the UV-IR relation ( |3.2| , p.3|) to rewrite /iexp (±z|0o) in terms 
of the IR parameters, the Hamiltonian can be made dimensionless h = H/M, depending 
only on the dimensionless volume / = ML, the coupling constant /? and r],{}. We compare 
the predictions on the spectrum, ground state energy etc. of the general two parameter 
BSG model to the truncated spectrum of this Hamiltonian. 



4.2 Finite size corrections from scattering theory 

Here we briefly recall the method to calculate the flnite size corrections for large volumes 
(/ ^ 1) from the knowledge of the bulk S'-matrices and boundary reflection factors. To 
simplify the presentation, let us consider a single scalar particle of mass m with reflection 
factors Ra {0) and RbiO) on the boundaries at x = and x = L respectively. Then the 
energy as a function of the volume can be obtained by solving the Bethe-Yang equation 

ml sinh ^ - i log i?„ (9) - i log Rb (9) = 2nl (4.3) 

for 9, where / is an integer (half integer) quantum number (corresponding to quantization 
of momentum in flnite volume). From the solution 9{L) of (|4.3| ) the energy with respect 
to the state with no particles is obtained as 

E (L) - E^'' (L) = m cosh 9 (L) . (4.4) 



Eqn. (|4.3| - |4.4| ) can also be used to give the {E{L),L) 'Bethe-Yang line' in a parametric 



form. When / = 0, eqn. (|4J3|) may have solutions corresponding to purely imaginary 
9, which may (in turn) correspond to boundary excited states obtained from the particle 
binding to one of the walls, cf. for details. 



4.3 Results 

In the TCSA for the general two parameter case the number of states with conformal 
energies below £^cut depends very sensitively on the coupling constant /3 (compactiflcation 
radius r), since the Hilbert space of the conformal free boson with Neumann boundary 
conditions is the direct sum of modules corresponding to the various momenta, which are 
integer multiple of 1/r. Therefore it is not a surprise that in the range > 3/2, where 
TCSA is expected to converge, there are so many states even for moderate -Ecut-s, that the 
time needed for diagonalizing H practically makes it impossible to proceed. 
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We overcome this difficulty partly by considering first only models on a "special line" 
in the parameter space described by 0o = or ^9 = 0. As pointed out in the models 
on this line admit the $ i— > — $ 'charge conjugation' symmetry as a result of the equality 
P"*" = P~ = P. As a consequence in these models there are two sectors, namely the 
even and the odd ones. It is straightforward to implement the projection onto the even 
and odd sectors in the conformal Hilbert spaces used in TCSA. This projection has two 
beneficial effects: on the one hand it effectively halves the number of states below -EcutQ, 
thus it drastically reduces the time needed to obtain the complete TCSA spectrum, and on 
the other the separate spectra of the even and odd sectors are less complex and therefore 
easier to study than the combined one. Furthermore, the spectrum of boundary states in 
the most general case depends only on t] [|], and so our considerations can be restricted to 
= without any loss of generality in this respect. 



4.3.1 Boundary energy 

First we investigate the ground state energy of these models to check the predictions of 
the BSG model. Since at one end of the strip we imposed ordinary Neumann boundary 
condition and switched on the boundary perturbation only at the other end, the ground 
state energy (in units of the soliton mass) for large enough L-s should depend on the 
dimensionless volume / = ML as 

where E{ri, {}) is the boundary energy, eqn. (p.4]) , and ?7Ar = f (1 + A) is the r] parameter of 
the Neumann boundary [Q. This prediction is compared to the TCSA data on Fig.s (|4.1[ - 
|4.2|) , where the dashed lines are given by eqn. ( |4.5| ). The agreement between the predictions 
and the data is so good that in the interval 5 < / < 15 the bulk energy constant and the 
sum of boundary energies can be measured with a reasonable accuracy. 

In our earlier paper P| when numerically investigating the ground state energy of the 
BSG model with Neumann boundary condition at both ends we made a conjecture that 

E{r]N,0) = -Ed{0) 

holds. {Ed is the boundary energy of the BSG model with Dirichlet boundary condition, 
eqn. ( |3.14|) ). Clearly the exact expressions eqn. ( |3.4| ) and eqn. ( |3.14|) do not satisfy this. 



but the violation of this relation is practically undetectable (using numerical methods) in 
the A range investigated in 0. 

We also checked the -d dependence of the boundary energy E{ri,{}), eqn. (0): the 
rapid growth in the number of states, caused by the absence of the two sectors, can be 
compensated by going to a sufficiently attractive value of A (A = 17) where TCSA is known 
to converge faster. In this case the choice i^cut = 13 resulted in 4147 conformal states and 

^In our numerical studies of these models E^ut varied between 15 and 18 and this resulted in 3 x lO'^ - 
5 X 10^ conformal states per sectors. 
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Figure 4.1: Ground state energy versus I in three BSG models with vq — y/An/P — 2 and 
1? = 0. 
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ure 4.2: Ground state energy versus / in four BSG models with 7] = O.Ttjn and 1^ = 0. 
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E{r]N, 0) + E{7], 13) (predicted) 


E{r]r,,0)+E{ri,^) (TCSA) 


5 


-0.22259 


-0.226959 


10 


-0.29012 


-0.29986 



Table 4.1: Boundary energies (in units of soliton mass) of two BSG models with A = 17 
and 1] = 0.7riN as measured from TCSA 



we could measure the sum of the two boundary energies fitting the volume dependence of 
the ground state energy by a straight line in the range 6 < / < 17, the results are collected 



in table 4.1 



To sum up, we showed that the prediction eqn. ( p.4| ) for the boundary energy of the 
general two parameter boundary sine-Gordon model is in perfect agreement with the TCSA 
data. This agreement indirectly confirms also the UV-IR relations, eqn. ( p^p.3D , since 
they were built into the TCSA program. The case of Dirichlet boundary conditions is 
investigated in the next section. 



4.3.2 Reflection factors and the spectrum of excited states 

We compare the reflection factors and the spectrum of excited states to the TCSA data 
in the case of models with 0o = (which is realized here as = 0). The bulk breathers 
naturally belong to one of the sectors, as the C parity of the n-th breather is (— 1)". 
However, since solitons and anti solitons can reflect into themselves as well as into their 
charge conjugate partners, solitonic one particle states (i.e. states, whose energy and 
momentum are related by i^^ = yP^-M^ where M is the soliton mass) are present in 
both sectors. 

To associate the various boundary bound states to the two sectors we have to determine 
the C parity of the poles z/„ and Wm in the soliton/antisoliton reflection factors. As in the 
even/odd sectors the reflection factors are given hj P±Q (where P = = P~ for ?9 = 0), 
the possible cancellation between the zeroes of Pq ± Qo and the poles of a{rj, u) have to be 
investigated. The outcome is that the poles at i'2k and W2k {k = 0, 1, 2, ..) appear in P + Q 
(i.e. the corresponding bound states are in the even sector), while the poles i'2k+i, W2k+i 
appear in P — Q (i.e. the corresponding bound states are in the odd sector). 

We analyzed the appearance of boundary bound states in the TCSA spectra of a number 
of BSG models. The results are illustrated on the example of a model when A = 7 and 
T] = 0.9 TjN. For these values of the parameters the sequence of Un-s and Wm-s in the 
physical strip is 

Uq > Wi > > W2 > 1^2 > W3 > 1^3 . (4.6) 

Therefore in the even sector we expect the following low lying bound states (i.e. ones with 
not more then three labels|): 

|0), |2), |0,2), |1,3), |0,1,1), |0,1,3), |1,2,3), |2,3,3), (4.7) 
^States having more labels are heavier thus they correspond to higher TCSA lines. 
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while in the odd sector 



|1), |3), |1,2), |2,3), |0,1,2), |1,2,2). (4.8) 

Since at one end of the strip the unperturbed Neumann boundary condition is imposed, 
the corresponding bound states are also expected to appear in the TCSA spectrum. As 
described in for r] = rj^ the z/„-s and the coincide and the bound states can be 

labeled by an increasing sequence of positive integers |ni, n^) n with Uk < A/2. Therefore 
in the even sector there should be TCSA lines corresponding to the 

|2)^, |1,3)^, |1,2,3)^, (4.9) 

'Neumann bound states', while in the odd one to 

|3);v, |l,2)iv, |2,3)^. (4.10) 

Finally there should be TCSA lines describing the situation when both boundaries are in 
excited states with no particle(s) moving between them, thus e.g. one expects a line in the 
even (odd) sector that corresponds to |0) ® |2)Ar (|0) ® |l)Ar)- 



We compare the predictions about these bound states to the TCSA data on Fig.(^) 
where the dimensionless energy levels above the ground state are plotted against /. On 
both plots the continuous lines are the interpolated TCSA data and the various symbols 
mark the data corresponding to the various boundary bound states and Bethe-Yang lines0. 

The two plots on Fig. (|4.3|) show in a convincing way that the low lying boundary 
states indeed appear as predicted by the bootstrap solution. (We show only those really 
low lying ones whose identification is beyond any doubt; the higher lying ones may be 
lost among the multitude of other TCSA lines). The reader's attention is called to two 
relevant points: first there is no TCSA line that would correspond to a |1, 1) bound state. 
The absence of this state is explained in the bootstrap solution by a Coleman-Thun 
diagram, that exists only if Wi > Ui. Second, both in the even and in the odd sectors, there 
is evidence for the existence of the lowest bound states with three labels. These states are 
predicted in the bootstrap solution by the absence of any Coleman-Thun diagrams when 
^ni > > hoMs. Thcsc two findings together give an indirect argument for the 
correctness of the boundary Coleman-Thun mechanism. This is most welcome, as the 
theoretical foundations of the boundary version of this mechanism are less solid than that 
of the bulk one. 

On the plots on Fig.(^) we also show in case of the lightest breathers B^, the 
excellent agreement between the TCSA data and the energy levels as predicted by the 
Bethe-Yang equations ( [4.3p.4|) , using either the ground state refiection factors (p.6| , [2?^ ) 
or the ones on the |0) excited boundary (|2.9| , |2.10|) . (In the latter case one has to take 



into account that now k is odd, b'^{ri,u) = 1, and the energy above the ground state also 
contains the energy of |0)). 

^°Some of the higher TCSA Hnes appear to have been broken, the apparent turning points are in fact 
level crossings with the other line not shown. This happens because our numerical routine, instead of 
giving the eigenvalues of the Hamiltonian in increasing order at each value of I, fixes their order at a 
particular small I and follows them - keeping their order - according to some criteria as I is changing to 
higher values. 
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Figure 4.3: TCSA data, boundary bound states and breather Bethe Yang lines in the BSG 
model with A = 7 and rj = 0.9 rj^. 
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A 


1^0 


-Ebuik (exact) 


£;buik (TCSA) 


-^boundary (cxact) 


-^'boundary (TCSA) 


31 





-0.01267857 


-0.01267(2) 


-0.0259997 


-0.026(17) 


31 


0.2 


-0.01267857 


-0.0126(14) 


0.1773231 


0.17(30) 


31 


0.495 


-0.01267857 


-0.012(25) 


1.009779 


0.97(75) 


17 





-0.02316291 


-0.0231(22) 


-0.0484739 


-0.049(06) 


17 


0.485 


-0.02316291 


-0.022(67) 


0.998483 


0.97(78) 


17 


0.5 


-0.02316291 


-0.022(69) 


1.048474 


1.02(84) 


7 


0.25 


-0.05706087 


-0.056(25) 


0.259213 


0.24(88) 


7 


0.48 


-0.05706087 


-0.055(62) 


1.054646 


1.03(23) 


41/8 


0.25 


-0.07911730 


-0.077(36) 


0.2464426 


0.23(14) 


41/8 


0.36 


-0.07911730 


-0.076(92) 


0.6381842 


0.61(72) 


41/8 


0.44 


-0.07911730 


-0.076(56) 


0.9513045 


0.92(52) 


7/2 





-0.1203937 


-0.118(22) 


-0.2957454 


-0.30(11) 


7/2 


0.3 


-0.1203937 


-0.114(69) 


0.4241742 


0.39(37) 


7/2 


0.42 


-0.1203937 


-0.11(34) 


0.9532802 


0.91(12) 


7/2 


0.5 


-0.1203937 


-0.11(18) 


1.295745 


1.23(60) 



Table 5.1: Boundary energy for Dirichlet boundary conditions: comparison to the TCSA 
data. The values for the boundary energy are for two identical boundary conditions at 
both ends of the strip. Energies are given in units of the soliton mass. 



5 TCSA: Dirichlet boundary conditions 

For Dirichlet boundary conditions, the formula ([4.1|) has to be changed: the terms con- 
taining boundary perturbations must be omitted, since there are no relevant boundary 
operators on a Dirichlet boundary. Furthermore, one must quantize the c = 1 free boson 
with Dirichlet boundary condition, which preserves boundary conformal invariance as well 
as the Neumann one. The Hilbert space is also changed, because there is a single vertex 
operator (the identity) living on the boundary, therefore it is essentially the same as the 

vacuum module of the chiral algebra (which in this case is the U{1) affine Lie algebra). In 
all numerical computations the truncation level was i^cut = 22, which corresponds to 4508 
vectors. 



5.1 Boundary energy 

Here we summarize the agreement between the formula (|3.14|) , first derived in ||12| and 
TCSA with Dirichlet boundary conditions. At both ends of the strip, identical boundary 
conditions are imposed. In this case, it is easier to vary the field value 0o: the interaction 
needs to be calculated for each given value of the sine-Gordon coupling parameter A only 
once. The agreement between the predicted values of the bulk and boundary energy and 
the TCSA vacuum energy levels is illustrated on Figure |5.1| , while numerical results are 



summarized in table 5.1 
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ML 



Figure 5.1: Comparing the predicted bulk and boundary energies to the TCSA data for 
Dirichlet boundary conditions. The dots are the TCSA data for various values of A and 
/o = ^5 while the lines are their predicted asymptotic behaviour for large volume. 
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Figure 5.2: Checking the reflection factors of -Bi, B2 and -B3 for A = 31 and ^ = 0.2. 
The dots show the one-particle energies predicted from the Bethe-Yang equations, while 
the continuous lines are the TCSA results. All energies are relative to the ground state 
and are in units of the soliton mass. 



5.2 Reflection factors 

Using the Bethe-Yang equations ( [4.3| , |4.4|) , we checked that the predictions for the energy 
levels from the ground state reflection factors are in excellent agreement with the TCSA 
data. Figure |5]^ is just an illustrative example; for all the other values of A and 00 in Table 



5.1 we had similar results. The deviations are partly due to truncation effects, but partly 



signal the fact that the Bethe-Yang equation only gives an approximate description of the 
finite size corrections. 

On Figure |573| , we illustrate how to obtain excited boundary states by analytic contin- 
uation of one-particle lines. 



5.3 Spectrum of boundary excited states 

We also performed an analysis of boundary excited states for Dirichlet boundary condi- 
tions. As there are two identical boundaries, the states come in doublets with symmet- 
ric/antisymmetric wave functions if the two boundaries are in a different state, and are 
singlets if the two boundaries are in the same state. There is also a selection rule due to 
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Figure 5.3: Boundary excited states at A = 17 and ^ = 0.485. The upper line is the 
1 = one-particle Bi line, including its continuation to imaginary rapidities, while the 
lower line is another portion of the imaginary rapidity continuation coming from another 
solution of the Bethe-Yang equations. The two lines together fit very well to the energy 
level doublet corresponding to the combination of a boundary in its ground state |) and the 
other in the excited state |0, 1), at least for sufficiently large values of the volume parameter 
I. 
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A 


Wi 

2-K 


El — Eq (predicted) 


El - Eo (TCSA) 


31 


0.495 


0.032428 


0.0323(62) 


17 


0.485 


0.099750 


0.0997(68) 


7 


0.48 


0.14349 


0.143(82) 


7 


0.45 


0.35711 


0.357(94) 


41/8 


0.44 


0.44675 


0.447(62) 


41/8 


0.36 


1.0035 


1.00(64) 


17/8 


0.4 


0.89148 


0.89(73) 



Table 5.2: Energy of the first boundary excited state as measured from TCSA 



a parity introduced by Mattsson and Dorey; namely, whenever the excited state of the 
left boundary has an even/odd number of indices, the right boundary also has even/odd 
number of indices, respectively. 

For the cases when (po = ^, the first excited state is expected to be degenerate with 
the ground state and this is indeed what we found within numerical precision. For the 
other cases, the energies of the first excited state are summarized in table This state 
corresponds to both boundaries being in the same excited state, so it must be a singlet 
and its energy with respect to the ground state (in infinite volume) is predicted to equal 

-.--»--»4("-i)— -(^^-^) 

We can measure this energy difference using the TCSA data. The results are illustrated 



in table 5.2 



For higher excited states one can introduce the notion of level. For a state labeled 
as \ni, . . . ,nk) it can be defined as the sum of the integers labels ^n^. It turns out 
that the energies are more or less hierarchically ordered and increase with the level. We 
considered excited states up to and including level 2 (the first excited state is at level 0) and 
found excellent agreement with the predicted spectrum apart from cases when the TCSA 
spectrum was too dense to come up with a meaningful identification of the TCSA data 
points with individual states. We also fitted them with analytic continuation of breather 
lines where this was possible, which also agreed very well with the TCSA data (see e.g. 
figure 



6 Conclusions 

In this paper we described an extensive verification of some results on boundary sine- 
Gordon theory, comparing numerical TCSA calculations to predictions concerning the 
spectrum, scattering amplitudes, boundary energy and the identification of Lagrangian and 
bootstrap parameters of the theory. We found an excellent agreement and confirmed the 
general picture that was formed of boundary sine-Gordon theory in the previous literature. 
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The main open problems are the calculation of off-shell quantities (e.g. correlation 
functions) and exact finite size spectra. While correlation functions in general present a 
very hard problem even in theories without boundaries, in integrable theories significant 
progress was made using form factors. One-point functions of bulk operators have already 
been computed using form factor expansions in some simple boundary quantum field the- 
ories p5| and one could hope to extend these results further. In addition, the vacuum 
expectation values of boundary operators in sine-Gordon theory are also known exactly 
[T^ . It would be interesting to make further progress in this direction. 

Concerning finite size spectra, there is already a version of the so-called nonlinear 
integral equation for the vacuum (Casimir) energy with Dirichlet boundary conditions 
IT^ , but it is not yet clear how to extend it to describe excited states and more general 
boundary conditions, which also seems to be a fascinating problem. 
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